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ABSTRACT 

We present Spitzer Space Telescope photometry of six short-period polars, EF 
Eri, V347 Pav, VV Pup, V834 Cen, GG Leo, and MR Ser. We have combined 
the Spitzer Infrared Array Camera (3.6-8.0 /xm) data with the 2MASS JHK S 
photometry to construct the spectral energy distributions of these systems from 
the near- to mid-IR (1.235-8 /mi). We find that five out of the six polars have 
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flux densities in the mid-IR that are substantially in excess of the values expected 
from the stellar components alone. We have modeled the observed SEDs with 
a combination of contributions from the white dwarf, secondary star, and either 
cyclotron emission or a cool, circumbinary dust disk to fill in the long- wavelength 
excess. We find that a circumbinary dust disk is the most likely cause of the 8fim 
excess in all cases, but we have been unable to rule out the specific (but unlikely) 
case of completely optically thin cyclotron emission as the source of the observed 
8/zm flux density. While both model components can generate enough flux at 
8/iin, neither dust nor cyclotron emission alone can match the excess above the 
stellar components at all wavelengths. A model combining both cyclotron and 
dust contributions, possibly with some accretion-generated flux in the near-IR, 
is probably required, but our observed SEDs are not sufficiently well-sampled 
to constrain such a complicated model. If the 8 /im flux density is caused by 
the presence of a circumbinary dust disk, then our estimates of the masses of 
these disks are many orders of magnitude below the mass required to affect CV 
evolution. 

Subject headings: cataclysmic variables (V834 Cen, EF Eri, GG Leo, V347 Pav, 
VV Pup, MR Ser) — stars: magnetic fields — stars: low-mass — infrared: stars 



1. Introduction 

Cataclysmic variables (CVs) are semi-detached binary systems consisting of a white 
dwarf (WD) primary star that accretes material from a low mass secondary star, with typical 
orbital periods of P or b < 1 day. Due to its large specific orbital angular momentum, the mass 
lost from the secondary star doesn't fall directly onto the WD. In most CVs, it instead settles 
into a disk around the primary star before losing enough angular momentum through viscous 
processes to finally accrete onto the WD. The release of gravitational potential energy in 
the disk causes it to be, by far, the brightest component of a CV over a wide range of 
wavelengths. On the other hand, polars (or AM Her stars, after the prototype) are CVs 
containing a WD with a strong (B > 10 MG) magnetic field. The matter lost from the 
secondary star is captured by the field lines and funnelled onto the magnetic pole(s) of the 
WD. The infalling matter is shocked as it nears the surface of the WD, radiating over most 
of the electromagnetic spectrum from X- rays to i n frared . Observational characteristics and 



types of CVs are extensively reviewed by IWarnerl (119951 ) . 



Polars are highly variable, both at a low level (usually ~l-2 mag) over their orbital 
periods as the geometry of the system changes along the line of sight, and on much longer 
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timescales with much greater changes in luminosity (up to ~4 mag). Until recently, polars 
were generally thought to hover around a maximum brightness, occasionally switching into 
relatively short low states, with no preferred minimum brightness for a given system. The 
cause of these low states is not fully understood, but is possibly due to the presence of 
star spots on the surface of the secondary star passing across the L I point and suppressing 
mass transfer into the Roche lobe of the WD (IHessman et al.l 120001 ). The resultant lowered 
accretion rate leads to a period of lower luminosity until the star spot has moved away from 
the LI point and normal mass transfer resumes. More recently, however, a number of polar s 
have been seen to hover in low states (IGerke &: Howell! 120061 ; lAraujo-Betancor et al.l 120051 ) . 
suggesting that the true cause of the accretion state changes in these CVs is more complex 
than previously believed. 

The secular evolution of CVs is driven by the loss of orbital angular momentum, which 
causes the separation of the two component stars (and, hence, their orbital period) to contin- 
ually decrease. While the secondary star radius decreases with mass-loss, this simultaneous 
draining of angular momentum from the system keeps it in contact with its Roche lobe, 
and enables the ongoing transfer of matter to the WD. Observationally, this results in CVs 
with shorter orbital periods having secondary stars with correspondingly later spectral types 
(i.e., lower masses). Eventually this ongoing mass-loss from the secondary star causes it 
to drop below the hydrogen-burning limit, and the secondary star becomes a degenerate 
brown-dwarf-like object. This change in the internal structure of the secondary star causes 
its radius to increase in response to further mass-loss. This, in turn, results in a reversal 
of the trend of orbital period change, leading to a gradual increase in both P or b and binary 
separation. This is bel ieved to be responsible for the observed minimum orbital period of 
CVs at P or b ~ 80 min (IPaczynski fc Sienkiewicall98ll ; iRappaport et al.lll982l ; iHowell et al. 
200lh . 



Population synthesis studies predict that the second ary stars in 
should have ev o lved p ast the hydrogen-burning limit (e.g 



70% of all CVs 



Kolb 


1993: 


Howell et al. 


1997) 



Littlefair et al.l (120031 ) reviewed all relevant observational data and analyses in the literature 
and concluded that while there is plausible indirect evidence for the existence of brown dwarf 
secondary stars in 39 systems, none of these had a reliable secondary star mass estimate or 
accurate enough spectral type to be certain. More importantly, they concluded that even 
if all of the potential candidates do indeed contain brown dwarf secondary stars, this only 
amounts to about 10% of all CVs - a huge deficit compared to the predicted 70% of CVs. 
Consequently, it has become very important to ascertain the number of systems, if any, 
containing brown dwarf secondary stars, in order to place limits on our theories of CV 
population age and evolution. 
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Here we present Spitzer Space Telescope photometry of six polars that were identified as 



havin g very short orbital periods (P OT h ^$ 90 min) and/or very red 2MA SS colors (jHoard et al. 



20021 ). This, along with other published evidence (e.g., compiled by iLittlefair et al.l 120031 ) . 
suggests that they are good candidates for containing brown dwarf secondary stars. Our 
Spitzer observations were intended to enable us to unambiguously classify the spectral type 
of the secondary sta r by isolating its contribution to the mid-IR flux. However, as shown in 



Howell et al.l (l2006al . henceforth, Paper I), even in the mid-IR the signature of the low mass 
secondary star in these polars is contaminated by excess emission beyond the level expected 
from the stellar components. In this paper, we explore in detail the possible origins of the 
mid-IR emission in magnetic CVs. 

Preliminary results for four of our targets (EF Eridani, GG Leonis, V347 Pavonis, 
and RX JO 154. 0-5947) were presented in Paper I. The first three are presented again here 
with more extensive modeling for completeness of our entire sample of polars, but RX 
J0154. 0—5947 is poorly studied and too little is known about its system parameters to 
warrant the more detailed modeling in this paper (in addition, it is the faintest of our polars 
- see Paper I). In addition, we present the first results from Spitzer mid-IR observations of 
the polars V834 Centauri, VV Puppis, and MR Serpentis. 



2. Observations and Data Reduction 



All of our m id-IR observations were obtained with the Infrared A rray Camera (IRAC; 
Fazio et al.ll2004l ) on the Spitzer Space Telescope (jWerner et al.l 120041 ) . A log of the obser- 
vations can be found in Table [TJ The images were reduced and flux-calibrated with the S12 
version of the Spitzer IRAC pipeline, and downloaded from the archive as Basic Calibrated 
Data (BCD). We obtained five BCD frames per object per wavelength channel, dithered 
with the medium-scale gaussian pattern. The IRAC channels have central wavelengths of 
3.6, 4.5, 5.8, and 8.0 fim for channels 1-4, respectively. The BCD images were corrected 
for array-location-dependency using the correction frames provided by the Spitzer Science 
Center. The frames were then combined with dual-outlier rejection to produce a mosaicked 
image for each channel using the Spitzer Science Center mosaicking software, MOPEX0. 

We then performed IRAF^I aperture photometry on the mosaicked fields using an aper- 



1 i.e., Mosaicking and Point-source Extraction; see Makovoz fc Marleau ( 20051 ) and the IRAC Data Hand- 
book for a comprehensive discussion of the MOPEX package and the dual-outlier rejection algorithm. 

2 The Image Reduction and Analysis Facility is maintained and distributed by the National Optical 
Astronomy Observatory. 
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ture size of 2.5 pixels. These results were subsequently corrected to an aperture size of 10 
pixels using aperture corrections derived from 46 isolated sources in our fields. The values 
of the aperture correction factors are 1.190, 1.199, 1.211, and 1.295 for channels 1-4, respec- 
tively. From experience we found that the point response function (PRF) for the IRAC arrays 
is not well-modeled and we obtained more consistent results using aperture photometry than 
via PRF fitting. Sky subtraction was tested in three different ways: (1) using a sky annulus 
around the target; (2) using a sky region a short distance from the target (only for the more 
crowded fields of EF Eri and VV Pup); and, (3) subtracting a median sky from the whole 
mosaic before performing the photometry. The three methods were found to be consistent 
with each other. We extended the wavelength coverage for our dat a set into the near-IR 
using archival data from the 2MASS All Sky Point Source Catalogue (ISkrutskie et al.ll2006l ). 
The 2MASS J, H, and K s magnitudes were converted to flux density in mJy using the 
absolute calibration zero points given in I Cohen et al.l (120031 ). All of the target photometry 
can be found in Tableland the resultant spectral energy distributions (SEDs) are shown in 
Figure [H 



2.1. Error Analysis 



The only pu blished attempts at rigorous error analysis for IRAC photometry are by 
Lacy et al.l (120051 ). who derived absolu te errors of 7%, 12%, 11%, and 12% for channels 1-4, 



respectively, and by lReach et a 



(2005), who carried out an absolute calibration of the IRAC 



instrument. The analysis by lLacy et al.l (120051 ) was calibrated relative to corresponding 
2MASS flux densities in an effort to obtain an absolute uncertainty for the IRAC colors. 
Those results therefore incorporate the uncertainties in both the 2MASS photometry and 
the IRAC colors for their objects. As the absolute calibration for IRAC is likely to change 
as more data are obtained, we performed our own error analysis to determine the relative 
errors for our data, which can then be adjusted in the future to incorporate the absolute 
errors as they become better known. 

We performed aperture photometry on a number of isolated stars in each field (including 
fields for targets other than our polars), covering a wide range of flux densities. This was 
carried out for both the mosaicked images and the five individual BCD frames that were 
used to create each mosaic. The standard deviation of the five individual BCD flux density 
measurements of each test star was then plotted against the final flux density of the star 
measured from the mosaic (see Figure [2]), in the expectation that the standard deviation 
would decrease with increasing brightness. We find that there is indeed a correlation between 
the standard deviation and the final flux density, but there is considerable scatter around the 
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best polynomial fit to the data (solid lines in Figure [2]). To allow for this we subtracted the 
polynomial fit from the data, calculated the standard deviation of the residuals, and added 
this to the best fit at each flux density to derive the total error estimate as a function of flux 
density for each IRAC channel (dashed lines in Figure [2]). We adopt these uncertainties for 
all future IRAC flux density measurements. 



3. Infrared Spectral Energy Distribution Models 

In Paper I, we briefly discussed the likely contributions to the IR SEDs of polars from a 
number of expected system components. In this paper, we present detailed physical models 
to compare with our observed SEDs. The available model components are described below, 
and we discuss the modeling results in §U 



3.1. Initial Values of Model Parameters 



Our basic model includes several components (see below), each specified by a number 
of parameters (e.g., size, temperature, mass, etc.). We have held constant the values of 
as many as possible of these parameters using information obtained from various literature 
sources for each target (see Table [3]). Known relationships between physical parameters were 
used to infer the values of additional parameters from an observationally determined one; for 
example, given a WD mass based on observational da ta, we used the mass-ra dius relation 
for a magnetic, carbon-core WD with 7 = 0.1 from ISuh fc Mathews! (120001 ) to calculate 
the corresponding WD radius. Radii and tempe ratur e s corr espo nding to a g i ven se condary 
star mass and spectral type were obtained from I Cox! (120001 ) and iPont et al.l (120051 ). In the 
absence of measured physical parameters in the literature, sub-stellar (i.e., L type) secondary 
stars were assumed to have M 2 = O.O8M , R 2 = O.IR ^ and T ? . =1500-2000 K (with higher 



temperature corresponding to earlier L spectral type; iKirkpatrick et al.l 119991 ). In addition 
parameters relating to the geometry of the CV (e.g., binary separatio n, Roche su r faces 



etc.) were calculated via the well-known Roche potential relations (e.g., lKopalHl959l . [1969 



19701 ) that depend primarily on component masses and the orbital period, and can be solved 
numerically by utilizing an iterative approach (see Figure [3]). 

The masses and radii of the secondary stars in our models are, to first approximation, 
consistent with filling their Roche lobes; however, whether or not the Roche lobe is filled is 
extremely sensitive to small changes in the mass or radius of both the WD and secondary 
star. We expect that the masses and radii of the component stars, which we have obtained 
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from the literature or calculated as described above, are not exactly known. Consequently, 
we have generally not been overly concerned with ensuring that the model secondary stars 
precisely fill their calculated Roche lobes for a given CV system geometry. In a couple 
of cases, the secondary star spectral type required in the model that best reproduces the 
observed SED is not consistent with the mass estimates from the literature. In these cases, 
we suggest that either the mass estimates are inaccurate or we are encountering the well- 
known iss ue of the wide range o f observed masses for each spectral type of CV secondary 
star (e.g.. lSmith & DhillorJ[l998h . 



3.2. The White Dwarf 

The contribution of the WD primary star is simulated using a simple blackbody curve 
corresponding to the temperature of the WD, scaled to the WD size and distance to the CV. 
Although, in detail, the SEDs of WDs are not identical to those of blackbodies, the falling 
Rayleigh- Jeans tail of the blackbody curve in the IR is a reasonable approximation to the 
long wavelength SED of a WD. In any case, the largest WD contribution is typically only a 
few per cent of the total observed flux density at the shortest IR wavelength (i.e., 2MASS J 
band), and falls off rapidly to longer wavelengths. 



3.3. The Secondary Star 



The contribution of the secondary star is represented by choosing an empirical infrared 
SED template composed of JHK S (2MASS) and 3.6, 4.5 5.8 and 8.0 /mi (IRAC) obser- 
vations of a number of M, L, and T dwarfs (IPatten et al.ll2006l ). The Patten et al. sample 
includes stars of the following spectral types: M3.5, M5, M5.5, M6.5, M7, M8, M9, L0, LI, 
L2, L3, L4.5, L6, L8, TO. 5, T4.5, T5, T7, and T8. We caution that the sample includes only 
field stars, rather than stars in binaries that are subject to irradiation and mass-transfer 
effects, and so may not exactly reliably reproduce the SED of a CV secondary star. In 
comparisons made, however, it has generally been show n that CV secondary s tars do not 
deviate significantly from their isolated counterparts (e.g. lSmith fc Dhillonlll998l ) . and there 
is a long history of using field dwarfs to approximate the secondary stars in CVs. We also 
point out that no models of CV secondary stars exist at IRAC wavelengths, and therefore 
we use these templates as a first approximation. We compiled the distances to each of the 
template stars from various literature sources in order to transform all of their observed flux 
densities to a standard distance of 10 pc. Almost all of the template stars have distances 
measured via trigonometric parallax. The most distant template star is located at ~ 27 pc, 
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but 60% of them are closer than 10 pc. 

We then performed a linear interpolation between adjacent observed spectral types in the 
M, L, and T groups, separately at each of the seven wavelengths, in order to "fill-in" the SEDs 
of the spectral sequence in steps of 0.5 in spectral type. The observed data for spectral types 
of M5.5, L2, and T5 from the Patten et al. sample were not used in the interpolation process 
because they showed anomolous flux densities at all wavelengths compared to the relatively 
smooth transitions between adjacent spectral types seen for the remaining stars. Possible 
explanations for this include errors in the published distances or spectral type classifications 
for these stars, or they may posssibly just be anomolous examples of their spectral type. We 
used interpolated SEDs to replace the observed SEDs for the three anomolous stars. 



In all cases, iPatten et al.l (120061 ) estimate the uncertainties in the template flux densities 
as less than 2%. Transformation to the standard distance will have increased this uncertainty 
by a few per cent, up to ~10-15% in the worst case. All of our final model flux densities 
will be transformed to d — 10 pc using the distances to each of the polars (either as a fixed, 
known value or as a free model parameter - see §3.ip in order to utilize these empirical low 
mass stellar template SEDs as a model component. 



3.4. Cyclotron Emission 



As described in Paper I, we expect from a theoretical basis that cyclotron emission 
should be relatively unimportant as a contributor to the mid-IR SEDs of our polars, especially 
at the longest wavelengths. There is also empirical support for this claim. For example, 
Dhillonl (119981 ) found "no evidence of cyclotron humps" in the .fT-band IR spectrum of MR 
Ser. The "worst case" scenario occurs for EF Eri, which has the lowest magnetic field 
strength of our samp le of polars, with the cyclotron fundamental occuring at A ~ 8 /mi. Yet, 
Ferrario et al.l (119961 ) found only very weak (amplitude < 1 mJy) cyclotron features at the 
red end of their 1-2.4 iim spectr um of EF Eri, a result that is also supported by the i^-band 
spectra of lHarrison et al.l (120041 ). who found only a very weak cyclotron hump longward of 
2 /zm, with stronger features only at shorter wavelengths. 

In order to include a model component representing cyclotron emission, we first cal- 
culated a representative profile as a function of wavelength. The shape of the cyclotron 
emission profile for a given field strength (which is a sum over harmonics) initially rises 
as v 2 before falling off rapidly with increasing freque ncy after the transition from optically 



thick to optically thin at some harmonic 1 < m < 10 (ILamb fc Masterdll979l ; ICropperlll990 
Frank et al.ll2002l ). All of our polars except EF Eri have the cyclotron fundamental short- 
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ward of A ~ 4-6 /im, and we included the first 20 harmonics (which extends well into the 
optical regime in all cases). The corresponding profile of the cyclotron emission is strongly 
affected by a change in the transition harmonic from a value of m — 1, but the difference 
in the observed profiles for values of m > 2 is relatively small. In general, the wavelength 
of the peak of the cyclotron emission component shifts shortward in proportion to 1/m as 
m increases. So, at small (large) values of m, there is more cyclotron emission at longer 
(shorter) wavelengths in the model SED. When m > 3, we expect to see only the initial 
rising (toward shorter wavelength) power law part of the cyclotron profile in the mid-IR 
SEDs of our polars. We will investigate two cyclotron emission scenarios in our models, one 
in which the transition harmonic is at m = 1 or 2 (i.e., approximately completely optically 
thin emission) and another in which m > 3 (i.e., partially optically thick emission). 

Electrons approaching the WD experience a successively stronger magnetic field as their 
distance decreases. Thus, it is more appropriate to consider the cyclotron emission as a sum 
over emission components for a range of field strengths. In order to approximate this effect, 
we recalculated the cyclotron profile in 1 MG steps from 1 MG to the specified surface 
field strength. Each profile was scaled by the distance-dependent factor described in §3.1 
of Paper I (to account for the change in electron density and temperature with distance 
from the WD). All of the profiles were then added together to obtain the final cyclotron 
emission profile for a given field. For the same surface field strength (B), the "sum over field 
strengths" prescription generally results in more cyclotron flux at longer wavelengths than 
in the single-field case. It also produces a cyclotron "pseudo-continuum," which tends to 
smooth over individual cyclotron humps (we already do not expect to see a strong signature 
of cyclotron humps in our model SEDs due to the poor spectral resolution of our data, which 
is matched in the model calculations). 

The final parameter for the cyclotron component is a scale factor, which is multiplied 
with the model profile to match the observed flux densities. This model component is the 
only non-physical component, in the sense that the absolute flux level is not predicted solely 
from other physical parameters (like temperature, distance, etc.). However, we are able to 
strongly constrain the scale factor for the cyclotron component, because when it is combined 
with the WD and secondary star components, the total cannot exceed any of the observed 
flux densities. 



3.5. The Circumbinary Dust Disk 



The presence of a combina tion of optically thick and optically thin dust was first pro- 
posed by iBerriman et al.l (119851 ) as an explanation for the IR luminosity in CVs. Initial 
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evidence for the likely presence of dust in several polars observed with Spitzer was described 
by us in Paper I. We now approach this analysis in more detail by including a realistic 
estimate of dust emission as a component in our SED model. 

There are several possible origins for dust in CVs, including (1) remnant dust produced 
during the common envelope phase and ejected from the binary with some angular momen- 
tum, (2) dust lost from the accretion flow and ejected from the binary, (3) dust produced in 
the mass outflows during r epeated nova or dw arf nova outbursts that does not escape from the 
gravity of the binary (e.g.. ICiardi et al.ll2006l ). and/or (4) dust produce d over long time-scales 



in the mass outflow of a wind from the accretion disk and / or W D (seelTaam fe Spruitl 12001 



Spruit fc Taamll200ll . iTaam et al.ll2003l . iDubus et all 12004 and IWillems et al 



2005 



for fur- 



ther discussion of the possible origin of circumbinary dust in CVs). 



Dust that is located outside the Roche lobes of the stellar components in a CV should 



settle into a rela t ively cool, geometr ically thin disk (ISpruit fc Taam 



2001; Taam et al. 2003 



Belle et al.l 12004 ; IDubus et al.l 120041 ) . There is a growing body of observational evidence for 
circumbina ry dust disks in both post-AGB binary s ystems (IVan Winckelll2004l ; lDe Ruyter et al 
20051 . 120061 ) and pre- main-sequence binary stars (jHartmann et al.l 120051 ) . The detection of 
crystalline silicates in these circu mbinary dust disks has been interpreted as an indicator 



of the long lifetimes of the disks (jHartmann et al.ll2005l ; iDeroo et al 



20061 ). As there is no 



previous observational evidence of circumbinary dust disks around CVs, it is impossible to 
say whether they have a similar longevity. However, the assumed temperature and density 
of the circumbinary dis ks are approximatel y an order of magnitude too low for significant 
sputtering of the grains (IBurke fc Silkill974l ). and if we assume that they are similar to those 
observed around other binary systems, then the dust disks will require little or no replen- 
ishment. The total minimum mass of dust that must be ejected from the system in order 
to form a circumbinary disk is therefore relatively small (approximately the mass of a large 
asteroid). This assumes that the dust is formed in the inner system pri or to ejectiom most 



likely in the outer atmosphere or wind of the low-m ass secondary star (ILunine et al 



1989 



Allard et al.l l200ll : iLeggett et al.l 120021 ; iTsuiil |2002|) - the probability of the dust forming 
in-situ in the disk itself is extremely low, as the dust formation mechanism requires high 
gas densities and, th erefore, would require a high mass transfer rate from the inner binary 
system into the disk (IBerriman et al.lll985l ). 



Recently, iMuno fc Mauerhanl (120061 ) reported the discovery of a Spitzer /IRAC mid-IR 
excess around several low mass X-ray binaries (LMXBs), which they interpret as evidence 
for circumbinary dust disks with almost identical characteristics to the ones we model in 
this paper. Searches for dust dis ks around CVs h ave so far failed to uncover any direct 
evidence for their existence (e.g., iBelle et al.l 120041 ) . In most CVs the circumbinary dust 
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disk flux in the IR (if present) is likely to be overwhelmed by radiation from the cool outer 
regions of the accretion disk (or even the bright, early spectral type secondary star in long 
orbital period systems). In that regard, polars may be the only type of CV in which the 
flux contribution from cool dust can be detected. Despite the lack of direct observational 
evidence, the existence of circumbinary disks in CVs has been suggested as the explanation 
for both the high mass transfer rate in some systems, as well as the observed l ack of pile- up of 



systems with orbital per i ods near the period m inimum (ISpruit fc Taamll200ll ; iTaam fc Spruit 



200ll ; ITaam et al.l 120031 ; IWillems et al.l 120051 ) . However, as described in more detail in £0 
the estimated total masses of our model circumbinary disks are likely too small to influence 
CV evolution as described in the papers cited above. 

For inclusion as a component in our SED model, we will initially make the assumption 
that the circumbinary dust disk is optically thick, a s this presents a computationally simpler 
case, and even (as described in iDubus et al.l 12004 ) might be the most appropriate case for 
dust disks in CVs. For completeness, however, we will also consider the case of an optically 
thin dust disk. 



3.5.1. The Optically Thick Case 
Our optically thick, geometrically thin dust disk model for the polars is based on the 



descri ption of the expected properties of circumbinary disks in CVs given in IDubus et al. 



(120041 ) . We do not consider the case for a warped or flared dust disk as our small number of 
data points cannot support the extra parameters required to model anything more complex 
than the geometrically thin case. We utilize the analytic solution for an optically thick dust 
disk r a dial temperature pro file and corresponding integrated flux density derived in 



(2003), iBecklin et al.l (120051 ). and references therein. In particular, equation 1 from 



Jura 



Jura 



(120031 ) provides the temperature profile of the disk: 



2 \ ( R *\ r 
-disk ~ I — -?*, 



1/4 / p \ 3/4 

3tt J 

where T* and R* are the temperature and radius of the central irradiating star, respect ively. 



and T rHsk is the resultant temperature of the disk at radial distance r. Equation 2 from I Jura 



( 120031 ) provides the integrated flux density, F disk : 



r ut 



„ 2,7V COS % , 

F diA =— ^— / B u (T disk )rdr, (2) 
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where i and d are the system inclination and distance, respectively, r- m and r out are the inner 
and outer radii of the disk, and B v is the Planck function. In this formulation, the dust 
disk is assumed to be optically thick but geometrically thin (i.e., with negligible height), 
so the integrated flux density originates only from the visible surface of the disk (with no 
contribution from the negligibly tall disk rim). 

In principle, the optically thick dust model is not sensitive to the exact properties of 
the dust grains (size, composition, etc.), as long as the grains are assumed to uniformly 
(re)radiate as blackbodies. The exponent on the radius-dependence in the equation for the 
disk radial temperature profile, however, does have some implication for the assumed prop- 
erties of the dust grains. A large exponent leads to lower temperatures near the inner edge 
of the disk, which corresponds to large dust grains that cool efficiently Conversely, a small 
exponent leads to higher temperatures near the outer edge of the disk, which corresponds to 
small dust grains that do not cool efficiently. In practice, we have found that an exponent of 
3/4, as used in equation [TJ is the only value that produces viable results. Larger exponents 
produce disk SEDs that rise too shallowly at longer wavelengths and are, overall, too faint 
to match the observed mid-IR flux densities without arbitrarily increasing the temperature 
of the disk's inner edge far beyond what could be produced by the irradiating stars. Smaller 
exponents produce disk SEDs that rise too steeply and are, overall, too bright to reproduce 
the observed SEDs without arbitrarily increasing the distances to the CVs to many hundreds 
or thousands of pc. Departure from an exponent of 3/4 implies non-blackbody grains, which 
are further discussed in Section 3.5.2. 



The disks considered by lJural (120031 ) and iBecklin et al.l (120051 ) are around isolated WDs 



(G29-38 and GD362, respectively). There are a number of differences that must be considered 
when applying this model to a circumbinary disk in a CV system. First, the disk will be 
centered around the center-of-mass of the CV instead of the WD. In the case of a low mass 
(potentially sub-stellar) secondary star, the difference between the center of the WD and 
the center-of-mass of the CV is relatively small (< 10% of the binary separation). Second, 
the inner edge of the disk will be much further from the WD. In the WD disk models, 
the inner edge of the disk is typically located at a few to ~ 10 WD radii. However, in 
the case of a circumbinary disk, the inner edge is, by definition, located outside the Roche 
lobes of the stellar components. We assume an inner radius for the circumbinary disk equal 
to t he tidal truncation radius at 1.7 times the separation of the stellar components in the 



CV (iDubus et al.ll2004l ). For our short-orbital-period polars, this typically corresponds to 
r in > 100-Rwd- The outer edge of the disk is calculated to correspond to a temperature of 20 
K using equation [TJ This typically corresponds to ~0.5-2xl0 4 WD radii for our polars, but 
in reality, the radius of the outer disk has a negligible effect on the model flux densities at 
the IRAC wavelengths, and serves only as a conservative, consistent cut-off value. Third, we 
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must consider irradiation heating of the disk from not only the WD but also the secondary 
staJE For example, for a typical short-period CV composed of a WD with temperature 
T WD = 10, 000 K and radius -Rwd = O.O1_R , and a low mass secondary star with T 2 = 3000 
K (i.e., 0.3T WD ) and radius R 2 = 0.1R Q (i.e., 10.Rwd), the ratio of total luminosities (given by 
L oc R 2 T 4 ) is L 2 /Lwd = 0.8. Thus, the total luminosity of the secondary star is comparable 
to that of the WD. 

Unfortunately, the additional contribution from the secondary star makes our consider- 
ation of the heating of the circumbinary disk more complex, because the secondary star is 
much further from the CV's center-of-mass (i.e., the center of the circumbinary disk) than 
the WD, so is, at any instant, closer to the inner edge of the circumbinary disk. However, 
under the assumptions that the disk inner edge radius is large compared to the binary sepa- 
ration (see above) and that the orbital period is short (so, on a time-averaged basis, all parts 
of the disk "see" the secondary star equally), we will neglect the offset of the secondary star 
and assume that all disk irradiating sources are located at the center-of-mass. 

Using the typical CV parameters listed above with P or ^ = 90 min, and stellar component 
masses of Mwd — 0.6Af Q and M 2 = O.1M , the circumbinary disk will have an inner edge 
radius of i?; n = 100-Rwd- Assuming only heating from the WD, equation [1] predicts an 
inner edge temperature of T; n ?s 200 K. If only the secondary star is considered as the 
heating source, then we get T in « 350 K. We have considered two methods for estimating 
the temperature of the circumbinary disk's inner edge that is likely to result from heating by 
both stellar components in the CV. The first method is to simply add together the inner edge 
temperatures predicted by heating from each of the stellar components individually. This 
would give T in rs 550 K for the typical CV parameters described above. The second method is 
to calculate a single "equivalent" star with radius -R eq uiv that is equal to the quadrature sum 
of R-wd and R 2 , and temperature T equiv such that L equiv = Rl quiv T^ quiv is equal to the sum of 
the individual luminosities of the WD and secondary star. For the system parameters listed 
above, this gives T m m 450 K. Given the uncertainty in estimating the likely heating effect of 
both stars, and the fact that both of our suggested methods yield similar temperatures, we 
have used the circumbinary disk inner temperature as a free parameter in our model SED 
calculations. However, we have constrained T in to be within fa ±200 K of the value predicted 
for heating by a single "equivalent" star (i.e., the second method described above). 



3 Fortunately, since polars lack an accretion disk, we do not have to consider this additional potential 
irradiating source. In principle, the accretion spot(s) at the WD magnetic pole(s) could be considered as 
additional discrete sources of irradiation; however, unless the magnetic axis is aligned very close to the WD 
equator, the geometrically thin circumbinary disk will not "see" the geometrically small accretion spot(s). 
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3.5.2. The Optically Thin Case 



The basic structure of an optically thin circumbinary dust disk is identical to that in 
the optically thick case - it is geometrically thin and heated by irradiation from both stellar 
components, both of which are assumed to be located at the system center-of-mass. Our 
model is based on two main assumptions - first, that each dust particle (re)radiates as a 
blackbody and, second, that all of the particles are "seen" by the observer (i.e., there is no 
eclipsing of the disk by the stellar components and every one of the grains contributes to 
the total flux density). The former condition is not an unreasonable assumption since the 
inner disk radius is set at the tidal truncation ra dius, well aw ay f rom the stellar compo nents. 
We use the same disk temperature profile from lJural (120031 ) and iBecklin et al.l (120051 ) as we 
use in the optically thick case, with an outer disk radius corresponding to an outer disk 
temperature of T out = 20 K, also as used in the optically thick case. We assume that the 
disk has a constant volume density throughout, so the volume of material at any temperature 
(radius) in the disk increases with falling tempe ratur e (incr easing radius). We used a fixed 
grain radius of 10 /im, which was suggested by lJural (120031 ) as the minimum grain size for 
dust debris disks around main sequence stars, and a fixed dust grain density of 3 gem -3 . 
In fact, the exact values of the dust grain parameters that we use have little effect on the 
resultant model SED of the dust disk; for example, if we change the grain size (thereby 
reducing the available radiating surface area in the disk), then the required luminosity to 
match the observed SED can still be achieved by changing one of the other disk parameters, 
such as the volume density of the dust. 

In all cases, we kept the disk thickness and volume density as low as possible to ensure 
that the disk remained optically thin; however, this was very difficult to quantify. The optical 
depth of the dust is given by: 



T 



K s pS, 



(3) 



where k s is the mass absorption coefficient, p is the volume density of the dust, and S is 
the column length along the line of sight. The mass absorption coefficient is a function of 
dust composition, grain size and shape, dust temperature, and the frequency of the incident 
radiation. As a result, it is usually quoted in a frequency-dependent form, k v . It is highly 
temperature-sensitive, but extremely poorly studied for temperatures above ~300K, as all of 
the literature on this subject deals with the cool interstellar medium and there are no pub- 
lished values of k v at the temperatures found in our disks. Even the theoretically calculated 
values for the relatively we ll-studied low temperatures vary by a factor of 30, from 0.002 to 
0.10 cm 2 g -1 (lDraindll989l ). As a result, we have no certain way of testing the optical depth 
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of our disks, but we adopted a simple, generic relation in an attempt to quantify our values: 



k v = 0.1(z//10 1J Hzf cm 2 g 



(4) 



which is the general form given in iBeckwith et al.l (119901 ) . The value of (3 depends on the 
growth s hape of the grains - spheriodal grains will have a higher value of (3 than fractal 



ones^see IBeckwith et al.l ll990L fo r an in-depth discussion of these effects). We follow the 
example of IBeckwith et al.l (119901 ) and accommodate the uncertainty in the grain shape by 
adopting the general relation given above for k v with (3=1. In any case, the volume density, 
and therefore the optical depth, of the disk is dependent on the disk temperature profile, 
the grain size and mass density, and the inner and outer radii of the disk. Any calculation 
of the true optical depth of the disk requires better constraints on all of these parameters. 

When the presence of dust can explain the long wavelength end of a polar SED, we 
generally have been unable to differentiate between the optically thick and thin circumbinary 
dust disks - because of the number of available parameters used to describe each type of 
disk, we can produce qualitatively similar model SEDs using either type of disk. Thus, in 
the remainder of this paper we will generally only discuss the less complex optically thick 
case, with the exception that we will consider an estimate of the total dust mass from the 
optically thin case (see §5]). 



4. Modeling Results 

We emphasize that we have seven data points in each of our SEDs, only four of which 
(from IRAC) are in the mid-IR. Consequently, our modeling results are intended to illustrate 
general trends in polars as a class rather than provide a specific, detailed representation of 
each of our target polars. In general, we started each model by including best estimates 
for the WD and secondary star components, then attempted to fit the observed 8 /im flux 
density through the addition of either cyclotron or dust emission, without exceeding the 
observed flux density at any shorter wavelength. In reality, these systems could (and, likely, 
do) include contributions from both cyclotron and dust emission, but we have too few data 
points to constrain such a complex model in any meaningful way. However, as will be 
described in detail in the following sections, success at reproducing the observed 8 /im point 
with a particular model configuration while simultaneously not exceeding the observed flux 
densities at shorter wavelengths can allow us to rule out cyclotron emission (often) or dust 
(occasionally) as a dominant contributor at the long wavelength end of the observed SEDs. 

Because of the relatively small number of observational constraints compared to the 
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degrees of freedom in the models, optimization of model parameters was accomplished via 
manual testing. In general, the values of the model parameters can be varied within a range 
of about ±5-10% and will still produce a model SED of comparable quality to the "best" 
cases discussed below. For convenience in discussing our results, we divided the targets into 
three classes based on observed SED morphology. Class I contains EF Eri and V347 Pav, 
which have SEDs in which the IRAC portion is brighter than the 2MASS portion. Class II 
contains VV Pup, whose SED is most similar in shape to a falling Rayleigh- Jeans-like tail 
and may be dominated by the secondary star. Class III contains V834 Cen, GG Leo, and 
MR Ser, which have SEDs in which the 2MASS portion is bright and falls steeply toward 
long wavelengths but the IRAC portion is relatively flat. We will now discuss our modeling 
results for the polars in each of these classes. 



4.1. Class I - Bright IRAC 

The two polars in this class have oddly-shaped SEDs compared to the other four polars 
in our sample. Specifically, the IRAC portion of the SED, while still approximately flat at 
the longest wavelengths, is brighter than the 2MASS portion (see Figured]). This possibly 
stems from another potential complication to the interpretation and modeling of our polar 
SEDs; namely, that polars are highly variable objects. Their brightness changes dramatically 
on irregularly-spaced long-timescales (weeks to years) between high and low accretion states 
(by up to a few magnitudes), and over their orbital period (by as much as 2 mag). In order 
to best compare the non-simultaneous 2MASS and IRAC measurements for our targets, we 
attempted to determine both the accretion state and orbital phase of each of the polars at 
the times of both sets of measurements. 

We used flux measurements from the literature and long-term light curves from the 
American Association of Variable Star Observers (AAVSO) to check the accretion state of 
the polars. Of our six targets, we were only able to definitively confirm that one (EF Eri) 
was in the same accretion state (low) during both the 2MASS and Spitzer observations. We 
suspect that at least three of the systems (V834 Cen, GG Leo, and MR Ser; i.e., the Class III 
objects) were in high accretion states during the 2MASS observations because of the steepness 
of the 2MASS portion of their SEDs, which may be indicative of the falling Rayleigh- Jeans- 
like tail of accretion-generated luminosity. The 2MASS-IRAC SED transition in VV Pup 
(i.e., the Class II object) is very smooth, which suggests that the data were obtained during 
the same accretion state (regardless of which state). Light curves from the AAVSO show 
that VV Pup was in a high state for the 2MASS observations. There is no ground-based 
coverage at the time of the Spitzer observations, but from the long-term light curve, it is 
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obvious that the system spends most of its time in a high state, and the closest data points 
show no indication of the system dropping into a low state at that time. 

So far, there are no published orbital light curves of polars in the mid-IR, so we have no 
evidence that the flux levels at the IRAC wavelengths vary with the same phasing and/or 
amplitude as in the optical or 2MASS wavebands, but for completeness we consider this 
possibility We used published ephemerides to determine the orbital phasing of the 2MASS 
and IRAC data sets for our targets, and then used published optical or near-IR light curves 
to estimate whether the orbital brightness states of those data sets were comparable. Both 
EF Eri and V347 Pav were observed by 2MASS and Spitzer at significantly different orbital 
phases that correspond to different orbital brightness states (in both cases, the 2MASS data 
were obtained during a faint state, and the IRAC data during a bright state). The other four 
targets were observed by both 2MASS and Spitzer at similar orbital phases that correspond 
to the same orbital brightness state, and so the results at all wavelengths should be directly 
comparable. In the absence of mid-IR orbital light curves, we have assumed that the ratio of 
bright-to-faint flux density for the IRAC data is the same as observed in the optical/near-IR. 
In any case, there are enough model parameters to compensate for different bright/faint scale 
factors (and, of course, if the mid-IR is shown in the future to not vary over the orbital period 
like the optical/near-IR, then this contingency can be discarded!). We note that even after 
attempting a scaling correction, the SEDs for EF Eri and V347 Pav are morphologically 
similar to each other but different from the other classes we have defined (see §4.1.11 and 
33X21 for details). They show flat IRAC SEDs like Class III, but the 2MASS portions of 
their SEDs resemble Class II. 



4-1.1. EF Eridani 



We selected EF Eri for o ur IRAC campaign due, in large part, to its extremely red 
(J — K s > 2.5) 2MAS S colors (|Hoard et al.ll2002l ). which exceed even the reddest IR colors 



for mid-L dwarfs (Leggett et al. 



2002j). This C V was discovered as an X-ray source by 



Gurskv et al.1 (119781 ) and subsequently identified by lGriffiths et all (119791 ) and I Williams et al. 



( 1979 ) as a probable AM Her type. The magnetic nature of the system was confirmed by 
Tapial (119791 ) from measurements of the linear and circular polarization, whic h he used to 
derive the very short orbital period of P or b ~ 81 min, which was later refined by lPiirola et al. 
<jl987h . 



EF Eri entered an extended low state in 1997, only emerging back into a high state in 
March 2006, and only remaining bright for a few w eeks before dropping back in to the low 
state. The low state optical spectrum obtained by IWheatley & Ramsayl (119981 ) displayed 
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only weak line and cyclotron emission, indicat ing; that the low state m ass transfer rate is 
extremely small. Modeling of this spectrum by lBeuermann et al.1 (120001 ) set a limit of M9V 
or later for the secondary star spectral type, and a distance estimate of d = 100 ± 20 pc 
(consistent with the upp er limit of d = 119 pc estimated from near-IR observations by 
Young fc Schneider! 1 19811 ). iHowell fc Ciardil (120011 ) detected the secondary star in K-band 
spectra with a spectral type of either M6V or L4-5, although it is impossible to tell whether 
the star has dropped below the hydrogen burning limit based solely o n spectral type , as the 
temperatures may depend on irradiation and evolutionary history. lHarrison et al.l (120031 ) 
presented simultaneous multi-wavelength optical and infrared photometry in which they 
believed they detected a bro wn dwarf secondary star in the H and K bands, but follow-up 
phase-resolved spectroscopy (lHarrison et al.l 12004 ) showed that the variability in both bands 
is dominated by cyclotron emission throughout the orbital cycle, most likely originating at 
both of the WD magnetic poles. 



Recently, IHowell et al.l (j2006d ) determined an extremely low mass for the secondary 
star in EF Eri (M2 = 0.055M©), suggesting that an L spectral type classification is, indeed, 
appropriate. At a distance of d « 100 pc, the secondary star in EF Eri cannot be earlier 
than L5.0, or the combined flux densities of just the WD and secondary star will significantly 
exceed the observ e d J-b and point. Our model for EF Eri uses the WD parameters from 
Beuermann et al.l (120001 ) . the secondary star parameters from IHowell et al.l (j2006d j. and 
assumes a spectral type of L5.0 for the secondary star (see Table [3]). We then combined this 
model with an additional component representing cyclotron or dust emission to explore the 
limiting cases for the ability of either mechanism to explain the observed long wavelength 
SED of EF Eri (without exceeding the short wavelength end). 

Figure H] shows representative model SEDs for EF Eri. The model containing a circumbi- 
nary dust disk (Figure Hk, and Model 1 in Table 0J can easily reproduce the observed 8 jum 
point, but significantly underestimates the observed flux densities at the shorter wavelengths 
(except at the shortest wavelength, J-band). When the cyclotron transition harmonic is set 
to m — 1, the cyclotron emission is essentially indistinguishable from the circumbinary disk 
component, so, in the case of completely optically thin cyclotron emission, we cannot differ- 
entiate between cyclotron or dust producing the 8 /im emission. However, neither mechanism 
(on its own) can produce the high levels of emission at shorter wavelengths without greatly 
exceeding the 8 fim point. We found this to be true in all of our polars, so we do not mention 
it again in the discussion of our modeling results, and defer discussion of completely optically 
thin cyclotron emission to §3 

Figure Hb (see Model 2 in Table H|) shows the model SED having a cyclotron component 
with transition harmonic at m — 2, which best reproduces the observed long wavelength 



-19- 



end of the SED, but still underestimates both the short wavelength end and the 8 /im 
point (even when the model 5.8 fim flux density is pushed to the upper uncertainty limit 
of the observed point at that wavelength, as shown in the figure). Switching the cyclotron 
component to increasingly optically thick configurations (e.g., m = 3; see Figure Hfc and 
Model 3 in Table H]) recovers some of the short wavelength flux, but at the expense of 
substantially underestimating the IRAC portion of the SED. 

Up to now, we were using the sum-over-field-strengths prescription for the cyclotron 
emission (see §3.41) . As shown in Figure HH (Model 4 in Table HJ), the single field strength 
approach can approximately reproduce the bright observed H band point (at the expense 
of over-estimating the J-band point) because the cyclotron humps are not blended together 
into a "pseudo-continuum" as in the sum-over-field-strengths approach. It is possible that 
if the intensity scaling of profiles produced further from the WD (i.e., at smaller magnetic 
field strengths) falls off more steeply than estimated in Paper I, then the final profile might 
be dominated by the emission produced very near the WD surface and, hence, retain the 
sharply peaked shape seen in Figure HH. A better reproduction of the EF Eri SED might be 
achieved through a combination of models similar to Models 1 and 4 (i.e., the SEDs shown 
in Figures and -d). However, even in this case, there would be some "missing" flux at 3.6 
and 4.5 /im, suggesting the need for either model components with different spectral profiles 
or an additional emitting system component. Adding more model parameters/components 
is not warranted by the ability of these sparse data to constrain complex models. We note, 
however, that some experimentation demonstrated to us that a simple blackbody component 
with T m 1000 K produces an emission profile that peaks around 4.5 fim and could fill in 
the EF Eri model SED as required. 

As noted in §4.1[ the 2MASS and IRAC data for EF Eri were obtained during different 
orbital bright ness states (fain t and bright, respectively). The low-state i^-band photometric 



light curves of lHarrison et all (120031 ) show that bright /-Ffaint ~ 1-67. Assuming that this ratio 
holds into the mid-IR, we scaled the 2MASS data for EF Eri by a factor of 1.67 to match 
the short wavelength portion of the SED with the IRAC data, and repeated our modeling 
trials. We scaled the 2MASS data instead of the IRAC data in order to avoid altering the 
observed value of the 8 /im flux density (which is the main target of our modeling efforts). 
Modeling these scaled data yields substantially similar results to those described above for 
the unsealed data (e.g., Figure [5] and Models 5-6 in Table HJ). The model parameters must 
be changed slightly from the unsealed versions to account for the difference in the scaled 
flux density levels; notably, the secondary star spectral type can be as early as L3.5 without 
exceeding the scaled J-band point (at d = 105 pc). Regardless, the general trends seen in 
the unsealed models are preserved in the scaled models. 
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4-1.2. V347 Pavonis 



V347 Pav was discovered by iPounds et al.l (|1993l) in the ROSAT Wide Field Camera 
survey. It was identif i ed as a CV bvlO'Donoghue et al.l (119931 ) and confirmed as an AM Her 
type by iBailey et al.l (119951 ). lO'Donoghue et al.l (119931 ) showed that the optical light curve 
is complex, with variability over the orbital cycle of approximately 1 mag, and displaying 
substantial variations from one orbital cycle to the next. Their photometry an d spectroscopy 
both s uggested an orbital period of P or b ~ 90 min, which was later refined by iRamsay et al. 
(l2004f ) to P orb = 90.082219(52) min. 



Bailey et al.l (119951 ) fit the optical circular polarization and light intensity curves with 



two linearly extended cyclotron emission regions on the surface of the WD, with a model 
polar field strength of B = 25 MG and an orbital inclination of i = 60°. This was constrained 



further by IRamsay et al.l (119961 ). who used the orbital period -secondary star mass relation 
for a main sequence star to predict M 2 = 0.16MJJ which, combined with the fact that they 
had failed to observe any eclipses of the system, allowe d them to pu t an u pper limit on the 
inclination of i < 78°. Subsequent Stokes imaging by iPotter et al.l (120001 ) constrained the 
inclination more closely to 64° < i < 72°, with magnetic field strengths of B = 15 and 20 
MG for the upper and lower accretion poles, respectively. 

At the nominal literature distance of d = 171 pc, the secondary star can be no earlier 
than spectral type M6.0, or the combined secondary star plus WD (the latter contributes less 
than 8% of the observed SED at all wavelengths) exceeds the observed J and H points. Figure 
[6^ (Model 1 in Table [5]) shows the model result with the addition of a circumbinary dust 
disk. As with EF Eri, we can easily reproduce the 8 /im point, but the model underestimates 
the 3.6-5.8 /im points. Using a cyclotron component with m = 2 (Figure Eb and Model 2 
in Table [5]), the model SED has the closest match to the intermediate- wavelength data, but 
is no longer bright enough at 8 fim to match the observations. Figure [6b (Model 3 in Table 
[5]) shows the m = 3 case - larger values of m cause the cyclotron peak to continue shifting 
toward short wavelengths until the cyclotron component appears as a Rayleigh- Jeans-like tail 
(similar in shape to the long wavelength end of the secondary star SED). Figure EH (Model 
4 in Table [5]) shows an m = 2 cyclotron component using the single-field approach instead 
of the sum-over-field-strengths. As with EF Eri, this approach allows us to reproduce the 



4 We have use d this mas s in our models to predict a secondary star spectral type of M6.0 and radius 
of R 2 = 0.21i? Q (|Coxll2000h . However, we note that this radius is substantially larger than the calculated 
Roc he lobe in V 3 47 Pa v assuming Mwd = 1.00M Q ; the secondary star radius of i?2 = 0.16i? Q for this mass 
from lPont et~ZI |2005l ) is in better agreement with the Roche lobe size. In any case, the model SED results 
are not very sensitive to changes in the mass or radius of the stars. 
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sharp peak in the SED, which for V347 Pav is located in IRAC channel 1. However, also 
as with EF Eri, this approach does not reproduce the 8 /zm point, suggesting a combination 
of Models 1 and 4 (i.e., dust and single- field cyclotron) as the best approach to reproducing 
this SED. 

We now consider the possibility that the 2MASS and IRAC data were obtained at 
di fferent orbita l brigh tness states. A high-state if -band photometric light curve published 



in 



Bailey et al.l (119951 ) shows that V347 Pav is variable by a factor of ~ 4 over its orbital 
period. As with EF Eri, we assume that this behavior is present in the mid-IR also, and 
attempt to correct the 2MASS and IRAC data sets to acccount for orbital variability. To 
this end, we scaled the 2MASS flux densities by a factor of 4.0. At the nominal distance 
of d = 171 pc, the scaled SED can accomodate a secondary star as early as M3.5V, whose 
SED exactly matches the 3.6 and 4.5 /zm points. However, as not ed earlier in th is secti on, 



the M6.0V secondary star suggested by the mass estimate from iRamsay et al.l (119961 ) is 
already at the limit of being too large for the calculated Roche lobe size, and an M3.5V 
star simply will not fit into the Roche lobe of a CV with this orbital period. Consequently, 
we have retained the M6.0 V seco ndary star and adjusted the distance to the lower limit of 



the lAraujo-Betancor et al.l (120051 ) estimate, d = 133 pc. Figure [7^ shows the scaled data 
with a model containing a circumbinary dust disk component (Model 5 in Table [5]). There 
is considerable unmodeled flux density in the 2MASS portion of the SED, but even if we use 
an M3.5V secondary star, the model SED still somewhat underestimates these data points. 
At all values of the transition harmonic m > 2, the cyclotron component cannot match 
the 8 (Ltm point without exceeding the shorter wavelength IRAC points. Figure [7b (Model 
6 in Table [5]) shows a model in which the cyclotron component is partially optically thick 
(m = 2). The peak of the cyclotron emission component shifts toward longer wavelengths, 
leaving a deficit of modeled flux density at 8 /zm, but suggesting a path towards recovering 
the previously unmodeled flux at the 2MASS wavelengths. 



4.2. Class II — Dominant Secondary Star 



4.2.1. VV Puppis 



VV Pup was discovered during the early 20th century (Ivan Gentlll93ll ) as a fai nt (V 



14.5- 18; iDownes et al.l l200ll ). rapidly periodic (P or b = 100 .4 min) varia ble star (jWalker 



19651 ). Based on photometric and spectroscopic observations, iHerbigl (Il960l ) suggested it was 
a binary whose emission lines originated on the brighter compon ent. Subsequ ently, VV Pup 
was identified as the third known example of the AM Her class (jTapialll977l ). 
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VV Pup has the distinction among our sample of polars of having the observed SED that 
is most reminiscent of the falling Rayleigh- Jeans tail of a cool, main sequence star. Figure 
(Model 1 in Table [6]) shows a model SED that utilizes the r e cent tig ht constraints on the 
WD and secondary star parameters provided by iHowell et al.l (l2006bl ) with a circumbinary 
dust disk to provide the additional flux density required to match the observed 8 /im point. 
However, the model SED underestimates the observed flux densities from the H-band to 
4.5/zm (IRAC channel 2). This model requires a distance scaling of d = 100 pc, which is 
slightly more than la away f rom the nominal valu e of d = 129^ P c from the trigonometric 
parallax distance reported in lHowell et al.l (l2006bl ). 



For all values of the transition harmonic m > 2, the peak of a cyclotron emission 
component shifts to successively shorter wavelengths (e.g., between the K s and IRAC channel 
1 points for m = 2), which causes the total model SED to exceed all of the shorter wavelength 
points when the cyclotron emission is scaled to fit the 8 /zm point. By moving the polar back 
to the nominal distance of d = 129 pc and including a cyclotron component with m > 2 
(e.g., see Figure [Hb and Model 2 in Table [6]), we can mitigate this problem and produce very 
good model SEDs (with better agreement at 3.6 and 4.5/xm and only a small flux density 
deficit at 8/xm). However, this form of the model is not very discriminatory, as the shape of 
the cyclotron component SED for m > 2 is very similar to that of the secondary star, so this 
model is basically equivalent to including only the stellar components and simply scaling the 
distance up or down until the best match to the observed data is obtained. 



4.3. Class III - Bright 2MASS, Flat IRAC 

We suspect that the polars in Class III were in high accretion states during the 2MASS 
observations due to the shape of the bright 2MASS portion of their SEDs, which could be 
indicative of the falling Rayleigh- Jeans-like tail of accretion-generated luminosity. If present, 
then this (unmodeled) component likely contributes at 3.6 and 4.5 /j,m in addition to dom- 
inating the 2MASS portion of the SED. We have not attempted to model this additional 
component because the exact characteristics of the accretion-generated luminosity are not 
known in detail, and would depend on a large number of (largely unknown or poorly con- 
strained) additional parameters (e.g., mass transfer rate, plasma flow and magnetic field 
geometry, accretion spot size and shape, etc.). The additional complexity that would be 
introduced by adding another model component is not warranted for our sparsely sampled 
SEDs. 
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4-3.1. V834 Centaun 



Mason et al.l (119831 ) identified a 15th magnitude, blue, emission line (H I, He I, II) star 
(subsequently named V834 Cen) as the optical counterpart of the variable X-ray source 
H1409-45/E1405-451 (I Jensen et al.lll982l ). Their time-resolved optical photometry revealed 
a period of P or b = 101.52(2) minutes with an amplitude of ~ 1 mag. In addition, rapid 
variability on timescales of 1-3 seconds was obse rved; this was later shown to be quasiper i- 
odic in nature ; with a coherence time of ~ 1 min (jLarssonlll985l . Il992l ; llmamura et al.ll2000l ). 



Wright et al.l (119881 ) detected V834 Cen in the radio, at a peak level of 35 mJy at 8.4 GHz. 



The radio flux showed variability on timescales as short as 1 min. Based on these charac- 
teristics (co upl ed with its X-ray emission, a n d circular and linea r polarization observed by 



Tapialll982l and IVisvanathan fc Tuohyj|l983l ). iMason et al.l (119831 ) suggested that this object 



was an AM Her type CV. 



Recently, iMennickent et al.l (120041 ) have matched the secondary star features in near- 
infrared (1-2.5 /jm) spectra of V834 Cen with a spectral type of M8±0.5V. An M8V 
star is quite faint, so we ha ve used a minimum distance of d = 121 pc in our models 
(jAraujo-Betancor et al.l 120051 ). Figured^ (see Model 1 in Table [7]) shows the model SED 
with a circumbinary dust disk, which reproduces the observed 8 /im point but underestimates 
the flux densities at shorter wavelengths. The secondary star would contribute significantly 
to the total model SED only at much smaller distances (d ~ 50 pc) than estimated in the 
literature. Figure [9b (Model 2 in Table [7]) shows the model with the circumbinary disk 
replaced by cyclotron emission with m — 2. The cyclotron profile peak is centered around 
3.6-4.5 /im, and even when the cyclotron emission is scaled up so that all of the 3.6 /mi flux 
density is accounted for, it still underestimates the 8 /im point. Figure (Model 3 in Table 
|7J) shows a largely optically thick (m = 6) cyclotron component. In this case, the cyclotron 
component can account for some of the flux density in the bright 2MASS portion of the 
SED, but still cannot reproduce the 8 /im point. 



4-3.2. GG Leonis 



GG Leo was discovered b y the ROSAT All Sky Survey and subsequently identified as a 
polar by lBurwitz et al.l (119981 ) . They derived the orbital period P or b = 79.87946(7) min from 
X-ray dip timings and concluded that the shape of the light curve is consistent with accret ion 
at only one of the WD poles. The optical light curve presented in iBurwitz et al.l (Il998l ) is 
variable by about 1 mag over the orbital period, and their photopolarimetry confirmed the 
one-pole accretion model and allowed them to derive a magnetic field strength of B ~ 23 MG 
in the accretion region. A fit to the low-state spectrum and the cyclotron emission yielded 
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WD parameters of O.8M < M WD < 1.2M and T WD = 8000 K, although the autho rs note 
that temperatures up to T WD = 12, 000 K cannot be excluded. iRamsay et al.l (120041 ) give a 
revised mass estimate for the WD of Mwd = 1.13±0.03 M based on X-ray data. 



Burwitz et al.l (119981 ) estimated the secondary star parameters from the orbital period 
and the radial velocity profile, adopting a radius of R2 ~ 0.13i2 Q , a mass of M 2 ~ O.O9M , 
and a spectral type of M5V or later. (We have used a spectal type of M7.5V, which corre- 
sponds to these mass and radius values, in our models.) They then used the non-detection 
of spectral features expected from the secondary star to put a lower limit on the distance of 
d = 100 pc, while the upper distance limit of d — 300 pc is based on the fact that any larger 
distances would require unrealistically high accretion rates. Unless the distance to GG Leo 
is substantially larger than 300 pc, the secondary star spectral type cannot be earlier than 
M4V or the model flux at 3.6 /im will exceed the observed value. In addition, no secondary 
star template from the Patten et al. sample (even at early spectral types) can reproduce 
the observed 2MASS SED without greatly exceeding the observed IRAC data. At a median 
distance of d = 200 pc, the M7.5V secondary star SED underestimates the 3.6 and 4.5 /im 
points by about 50%, leaving room for the unmodeled accretion-luminosity component. 

As shown in Figure [TOa (Model 1 in Table [8]), a circumbinary dust disk can reproduce the 
8 /im point. A similar result can be obtained even if we assume the minimum distance (d = 
140 pc) at which the M7.5V secondary star exactly matches the 3.6 and 4.5 /im points. For 
all values of the transition harmonic m > 2, the peak of the cyclotron component is located 
at short wavelengths, making it impossible to reproduce the 8 /im point without greatly 
exceeding the 3.6-5.8 /im points. If we use only the single-field cyclotron prescription, and 
assume completely optically thick emission (m = 10), then the resultant model trends toward 
reproducing the bright 2MASS portion of the SED - suggesting that cyclotron emission, 
instead of or in addition to, accretion luminosity could also be responsible for the 2MASS 
data - but still cannot reproduce the 8 /im point (see Figure [TUb and Model 2 in Tabled]). 



4.3.3. MR Serpentis 



MR Ser was i nitially discovered in the Palomar-Green survey for UV excess objects 



(IGreen et al.lll982 ); it tu rned out to be the only magnetic CV discovered in the PG sur- 
vey. iLiebert et al.l (119821 ) performed the first detailed observational study of MR Ser, and 
measured an or bital period of - P nr h = 113.56 min from polarimetric observations. This was 
later refined by lSchwope et al.l (Il99ll ). who used time- resolved spectroscopic observations of 
a narrow emission line com ponent (believed to ori gi nate on the secondary s tar) t o measure 
P or b = 113.47 min (also see ISchwope et al.l Il993bh . IWickramasinghe et al.l ( I199ll ) observed 
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shallow cycl otron humps in their optical spectra, which varied with phase in amplitude and 
wavelength. lHarrison et al.l (120051 ) found cyclotron emission from the B = 26 MG WD mag- 
netic field contaminating the blue end of their i^-band spectrum, but the red end of the 
spectrum was not significantly affected by cyclotron emission. 

MR Ser dropped into a low state in 1985, at which time the optical s pectra showed wea k 
H-alph a emiss ion and absorption features of an M dwarf secondary star (jMukai et al.lll985l ). 
Szkodyl (119881 ) found that even in the low state, H and H e emission lines were st ill present, 
indicating that there was still some ongoing accretion. ISchwope et al.l (jl993bl ) identified 
features of the secondary star in long w avelength optical spectr a of MR Ser, corresponding 



to a s pectral type of M5-6V (also see iMukai &: Charles! Il986l ). However, lHarrison et al. 
( 120051 ) estimated a somewhat later secondary star spectral type, M8V, from an infrared 
spectrum. 

Of the three polars in our Class III, MR Ser and V834 Cen have the most similar 
observed SEDs. It is no surprise, then, that our modeling results for MR Ser are very similar 
to those for V834 Cen. We utilized a minimum distance of d = 134 pc (see Table [3]) in order 
to maximize the contribution from the faint secondary star. At this distance, the spectral 
type of the secondary star cannot be earlier than about M5.5V, or the model 3.6 /im point 
will be brighter than observed. Figure [TTa (Model 1 in Table [9]) shows the model for MR 
Ser containing an M7V secondary star (consistent with the mass estimate of M2 = O.1M 
- see Table [3]) and a circumbinary dust disk. As with the similar model for V834 Cen, this 
reproduces the observed 8 fim flux density, but underestimates the flux density at all shorter 
wavelengths. As also found for V834 Cen, at values of m > 2, the cyclotron profile peak 
shifts toward shorter wavelengths, making it impossible to reproduce the 8 /im point without 
exceeding one of the shorter wavelength IRAC points (e.g., see Figure [Tib and Model 2 in 
Table [9J . 

An interesting aspect of the observed IRAC SED for MR Ser is that the 4.5 and 5.8 /zm 
points are brighter than the 3.6 and 8.0 ^m points. This behavior is not seen in any of the 
other polars in our sample (see Figured]). Reproducing this shape is not possible using only 
circumbinary disk and/or cyclotron components, without invoking an additional source of 
flux that peaks between IRAC channels 2 and 3; for example a blackbody component with 
T fa 800 K (also see Appendix EJ). 
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4.4. Sample "Best" Models 

While keeping in mind our warning about the (in) ability of our data to constrain complex 
models, we also recognize that it is somewhat unsatisfying to just consider models that are 
designed to only reproduce the observed 8 [im point with little regard to the SED shape 
at shorter wavelengths (other than to avoid exceeding the observed flux densities). To that 
end, we will present a sample "best" model for each of our polars (see Table [TUI) . These 
models are not constrained to be unique solutions and should not be considered as definitive 
physical representations of the polars. They simply illustrate that it is possible to create 
reasonable reproductions of the observed IR SEDs for our polars, which, given additional 
data, could be developed into much more physically detailed models. However, the models 
shown in this section do not demand that the observed SEDs must be fit using only the 
model components utilized here. 

In all cases the "best" models involve combining the stellar components with cyclotron 
and circumbinary disk emission. Typically, the circumbinary disk emission dominates in 
producing the excess 8 fim emission, whereas the cyclotron emission dominates at shorter 
wavelengths. Figure [T2l shows sample "best" models for EF Eri (with unsealed 2MASS data), 
V347 Pav (unsealed), and VV Pup. In the case of EF Eri, we have optimised the model to 
best reproduce the IRAC SED, while neglecting the bright 2MASS H-band point (see section 
4.1.1 for an example of fitting the H-band point). Figure [T31 shows sample "best" models for 
V834 Cen, GG Leo, and MR Ser. In these cases, we have not attempted to reproduce the 
bright 2MASS portions of the SED. We can generally achieve good reproductions of the IRAC 
portions of the SEDs. Of course, there is likely contamination from the accretion-generate 
luminosity in IRAC channels 1 and 2, which, if modeled, would cause these "best" models to 
overestimate the flux density at those wavelengths. However, this can easily be accomodated 
by making small adjustments to the distance, circumbinary disk inner temperature, and/or 
cyclotron emission strength. Out of all six of our polars, the "best" model for VV Pup comes 
closest to exactly reproducing the observed SED (see §4. 2ft . 



5. Discussion and Conclusions 

We emphasize again that the results of our modeling should be considered less as specific, 
detailed pictures of each individual polar, and more as an illustration of the general mid-IR 
properties of the class of polars as a whole. With that in mind, we can describe the general 
trends that our modeling efforts have revealed: 

(1) There is excess flux density at 8 /mi above that expected from just the stellar 
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components in all of our target polars. Our capability to obtain sensitive IR observations 
has improved dramatically in recent years, especially following the launch of the Spitzer 
Space Telescope. Evidence for IR excess in isolated WDs is becoming increasingly common, 
and has been convincingly linked to the presence of circumstellar dust disks. Now, CVs and 
LMXBs have joined the pantheon of systems containing compact objects that show excess 
luminosity at long wavelengths. Among the CVs, polars (which lack the bright accretion 
disk of non-magnetic systems) may offer the best opportunity to detect the comparatively 
faint signature of circumbinary dust. The likelihood of detecting the dust emission is also 
improved by observing the CVs (polars or otherwise) with the shortest orbital periods (i.e., 
-Porb ^ 90 min), which contain the lowest mass (hence, least luminous) secondary stars. 

(2) In all circumbinary dust disk can reproduce the observed 8 fim flux den- 

sity. Equivalent results can be obtained using either optically thick or optically thin disks. 
Although we have only presented model SED results from the computationally simpler op- 
tically thick disks here, we have used equivalent models containing optically thin disks to 
provide an estimate of the masses of dust required to explain the observed mid-IR excesses 
in our polars. Typical values for the total dust mass in our optically thin circumbinary disk 
models are Md ust ~ 10 15 -10 17 g, or about 1O~ 19 -1O~ 17 M . This is only a very small frac- 
tion of the annual mass transfer budget in these short-orbital-period CVs, at typical rates 
of M < Kr n M yr" 1 Jflowell et al.lll995h . It is also smaller t han the 10 20 -10 22 g of dust 



produced during a superoutburst of a short-orbital-period CV (jCiardi et al.l 120061 1 . These 
circumbinary disks could be created (and, if necessary, replenished) from only a comparative 
trickle of outflowing matter over the course of millions or billions of years. 

In addition, the total dust mass in our mo del circumbinary d isks is many orders of 
magnitude lower tha n the ~ 10 29 g of ma t erial t hat lTaam et al.l (120031 ) require to significantly 



affect CV evolution. iMuno fc Mauerhanl (120061 ) also found a paucity of circumbinary material 
in their LMXB systems (Md ust ~ 10 22 g), as did lDubus et al.l (120041) from their circumbinar y 
disk model calculations for CVs (Md us t < 10 24 g). As noted by IMuno fc Mauerhanl (120061 ). 
even if the gas-to-dust mass ratio in the circumbinary disk is 100:1, then t h e tota l disk mass 
is still several orders of magnitude too small to satisfy the iTaam et al.l (120031 ) condition. 
If our estimates of the circumbinary disk masses are correct, then either they cannot be 
res ponsible for the ob served discrepancies in CV angular momentum loss rates as suggested 
by lTaam et al.l (120031 ). or the interaction between the circumbinary disk and the inner binary 
system is poorly understood. 

(3) On the other hand, a completely optically thin cyclotron emission model compo- 
nent typically has a calculated SED shape that is indistinguishable from, or very similar to, 
that of the circumbinary dust disk. Hence, we cannot exclude the specific case of cyclotron 
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emission with transition harmonic m = 1 as the source of the excess mid-IR flux density. 
However, when the model cyclotron emission is completely optically thin, it also does not 
contribute strongly at shorter wavelengths and the resultant model SEDs underestimate the 
flux density at short wavelengths. In addition, observational evidence for completely opti- 
cally thin cyclotron emission has not yet been found in any polar. Since the circumbinary 
disk emission always contributes most strongly at 8 /im (with little or no flux density con- 
triubution shortward of 3.6 /im), it seems more likely that cyclotron emission in these polars 
is not completely optically thin, and can account for some of the shorter wavelength flux 
density (e.g., see the "best" models discussed in §4.41) . 

(4) In all but one case, cyclotron emission with m > 2 is too faint at 8 /mi to reproduce 
the observed data when it is scaled so that it does not exceed any of the shorter wavelength 
IRAC points. The exception to this is VV Pup, which has the smallest mid-IR excess 
compared to the model stellar components. For this polar, optically thick cyclotron emission 
in combination with the stellar components can adequately reproduce the observed SED. 
However, this result for VV Pup does not really imply a preference for cyclotron emission 
since the optically thick cyclotron SED mimics the shape of the falling Rayleigh- Jeans tail 
of the secondary star SED. Hence, this model devolves to the equivalent of using only the 
secondary star model SED and scaling the distance to match the observed data. There are 
several possible resaons why VV Pup could contain the apparently smallest mid-IR excess 
of our sample of polars. It is probably not due to masking of the dust signature by a 
bright secondary star, since VV Pup does not have the longest orbital period of our sample 
(and, hence, does not have a correspondingly much earlier spectral type secondary star than 
the other polars - both MR Ser and V347 Pup have somewhat longer orbital periods and 
presumably similar secondary stars to VV Pup). If the IR excess is caused largely or solely 
by dust, then VV Pup possibly just implies that the amount or emissive characteristics of 
dust is not the same in all polars. The most likely explanation, however, is that VV Pup has 
the highest estimated inclination of our sample of polars, so the projected radiating area of 
the circumbinary disk is smallest. 

In summary, while a circumbinary dust disk can always explain the observed 8 /j,m 
excess in our target polars, only the specific case of completely optically thin cyclotron 
emission (which mimics the SED shape of a circumbinary disk but has not yet been observed 
in any polar) can also independently explain the mid-IR excess. However, even in this 
special case, the cyclotron emission can not be used to account for excess flux density at 
shorter wavelengths (where any feasible circumbinary disk emission is also negligible). In 
the presence of cyclotron emission that is at all optically thick (and, hence, peaks at shorter 
IR wavelengths), circumbinary dust disk emission would still have to be present in order to 
explain the 8 /im excess in our polars. If the mid-IR excess in our polars is produced, either 
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solely or in part, by dust in circumbinary disks, then the implied m ass of dust i s too s mall 



by many orders of magnitude to affect CV evolution as described in iTaam et al.l (120031 ) . 
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A. On the Possible Presence of a White Dwarf Circumstellar Dust Disk in 

Polars 

The SED shape produced by a single-temperature blackbody component can also be 
produced by an optically thick circumstellar dust disk arou nd the WD in a polar. Such a 



disk is exactly analogous to the WD debris disks described in Uural (120031 ) and iBecklin et al. 



( 120051 ). which were used as the basis for our consideration of the circumbinary disks. While 
the presence of circumstellar disks around isolated WDs is becoming increasingly "normal" 
as detailed IR observations of these objects are acquired, the existence of a WD circumstellar 
disk in a CV would be a rather novel occurence. Nonetheless, we will briefly explore the 
nature and plausibility of such a system component (although further exploration of its 
characteristics - and even its existence - will demand more detailed mid-IR observations 
than are presently available). 

The circumstellar disk is centered around the WD, and extends from the sublimation ra- 
dius close to the WD (i.e., the radius at which the disk temperature, according to equation [TJ 
is in the range 1000-2000 K, typically ~5-10-Rwd), out to the tidal tr uncation radius of 0.49 



times the binary separation («80% of the WD Roche lobe radius; e.g., iPapaloizou fc Pringle 



19771 ). The overall hotter temperature of this circumstellar disk can produce emission at the 
short IRAC wavelengths with an SED shape comparable to a single-temperature blackbody 
component, while its smaller radiating area (compared to a circumbinary disk) does not re- 
quire moving the CV to large distances in order to avoid exceeding the observed flux density 
levels. 
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It is important to clarify the nature of this circumstellar dust disk. It is more analogous 
to a passive planetary ring system (a la Saturn) than to a normal CV accretion disk. For 
example, it is assumed to be a uniform disk that is geometrically thin, with a typical thick- 
ness comparable to the size of the dust grains. Collisional interactions between dust grains 
are assumed to be relatively rare (compared to a viscous accretion disk). Consequently, we 
assume that there is no (or only negligible) accretion-generated luminosity in the circumstel- 
lar dust disk. Several potential sources of dust in a CV are listed in §3.51 We now add two 
more possibilities that apply particularly to dust in a circumstella r dis k around the WD: 



a tidally disrupted asteroid or comet, as described in I Jural (120031 ) and lBecklin et al.l (120 05). 



and (2) dust produced in the outer atmosphere or wind of the secondary star (ILunine et al. 



19891 : lAUard et al.ll200ll ; iLeeeett et al.ll2002t lTsuiill2002f ) and transferred into the Roche lobe 
of the WD through the LI point. 

In the polars, the WD magnetic field prevents the formation of an accretion disk because 
the plasma emerging from the LI point is captured onto magnetic field lines and funneled 
directly to the poles of the WD. This mechanism should not work efficiently to prevent the 
formation of a dust disk that is composed of un-ionized (or otherwise uncharged) material 
that will not interact strongly with the magnetic field. However, dust grains in space that 
are exposed to UV radiation will generate a flux of photo electrons and develop a positive 
surface charge that rapidly reaches an equilibrium value (jHoranyil Il996l ). The amount of 
equilibrium charge has some dependence on the grain composition, and the grain size deter- 
mines the timescale over which the equilibrium charge is reached (with larger grains reaching 
the equilibrium charge faster). We define the critical surface charge, Q C rit, on a dust grain 
such that, for larger charges, the dust will interact with the WD magnetic field in a polar 
and formation of a dust disk would be impeded in much the same way as for a normal CV 
accretion disk. For smaller charges, the dust and magnetic field will not interact appreciably, 
and the circumstellar dust disk can form even around the magnetic WD in a polar. 

By comparing the gravitational and magnetic forces acting on a dust grain emerging 
from LI and undergoing free-fall acceleration toward the WD, we can estimate the value of 
Qcrit- The gravitational force, F G , is given by 



GmM WD 



G 



(Al) 



where m is the mass of a dust grain, Mwd is the WD mass, r is the distance from the center 
of the WD to the dust grain, and G is the gravitation constant. The magnetic force acting 
on the dust grain, F m , is given by 
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where Q is the surface charge on a grain of dust, v is the free-fall velocity of the dust (itself 
a function of distance from the WD), and -Bwd is the WD surface magnetic field (magnetic 
flux density). After equating these two forces, we can solve for Q cri t to obtain 



Q 



'}GM WD ) 1/2 mR 3 s il 



crit 



(A3) 



where we have considered the maximum potential interaction between the WD magnetic 
field and the dust by solving for Q cr it as close as possible to the WD, at the dust sublimation 
radius, -R su b- 

We have calculated values of Q crit using the system parameters of MR Ser, and assuming 
-Rsub = 7i?wD, which corresponds to T su b = 1500 K (the midpoint of the possible range of 
sublimation temperatures). For spherical dust grains with uniform density of 3 g cm -3 and 
ra dii of 1, 10, an d 100 u rn, we obtain Q crit « 1.9 x 10 n e, where n — 2, 5, and 8, respectively 
Q. iKempf et al.l (120041 ) measured the charge on in-situ interplanetary dust grains in the 
solar system using the Cosmic Dust Analyzer on the Cassini spacecraft. They found values 

of ~8-34xl0 3 e for ~ 10 /xm, spheroidal grai ns. This is ~5-25 times smaller than our 

Qcrit for similarly sized grains. ISickafoose et al.l (120001 ) have experimentally investigated the 
photoelectric charging of dust grains, and find an equilibrium surface charge of ~ 4 x 10 4 e for 
graphite grains with radii of ~ 50 /zm in the presence of UV illumination with a minimum 
wavelength of ~ 2000A. Our v alue of Qc. r ;t for a grain o f this size is ~ 2 x 10 7 e, 500 times larger 
than the charge measured by ISickafoose et al.l (120001 ). We infer that the equilibrium charge 
that is likely to form on dust grains in WD circumstellar disks in our polars is substantially 
smaller than the charge required to produce significant interaction between the dust and the 
WD magnetic field at distances from the WD larger than the sublimation radius. Thus, the 
presence of a circumstellar dust disk around the WD in a polar, while possibly improbable, 
is not an impossible scenario. 

This calculation was motivated by the fact that, during our modeling process, we noticed 
that the unusual shape of the IRAC SED of MR Ser matches very closely to that expected 
from only a WD circumstellar disk. If the secondary star in MR Ser actually contributes very 



: We note in passing that our calculated values of Q C rit are several orders of magnitude smaller than the 
corresponding values, Qmax, at which the electrostatic tensile stress in the d ust grain interiors woul d be 
sufficient to fracture (i.e., destroy) the grains (e.g., following the calculation in Draine fc Salpeter 1979h . 
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little to the total SED (or at least contributes more or less equally in all IRAC channels), 
then the observed mid-IR SED can be reproduced very well through the addition of only 
a WD circumstellar disk with an inner temperature of 1300 K. However, in this case, the 
secondary star must be exceptionally faint in order to prevent the 3.6 /xm flux density from 
becoming too bright. Even an L type star produces too much flux at 3.6 jum at the distance 
required to match the flux density of the WD circumstellar disk (which is limited by the size 
of the WD Roche lobe) to the observed IRAC SED. A model with a T3.5 secondary star and 
a WD circumstellar disk with inner temperature of 1300 K produces an excellent fit to the 
IRAC SED, but only if the distance to the system is a mere 35 pc. This distance is a factor 
of y 4 smaller than the minim um distance quoted in the literature (which was calculated 
by lAraujo-Betancor et al.ll2005l from fitting the UV spectrum of the WD in MR Ser), not 
to mention requiring a much later spectral type for the secondary star than estimated in 
the literature. Consequently, we caution that, although this circumstellar-disk-dominated 
model is suggestive, inasmuch as it reproduces the observed IRAC SED so exactly, it is 
overly complex for the ability of our data to constrain it and is highly speculative. 
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Table 1. Log of observations with the Spitzer Space Telescope 



Polar 


AOR Key 


Observation Date 


Total Exposure Time 






(UT) 


(s) 


V834 Cen 


10191360 


16 Jul 2005 


20 


EF Eri 


10185984 


17 Jan 2005 


300 


GG Leo 


13532160 


06 May 2005 


120 


V347 Pav 


13531648 


19 Aug 2005 


300 


VV Pup 


10191872 


26 Nov 2004 


120 


MR Ser 


10188288 


27 Mar 2005 


20 



Table 2. Measured 2 MASS and IRAC flux densities (mJy) 



Polar 


J 

(1.2 fan) 


H 

(1.7 /im) 


K s 
(2.2 fim) 


IRAC-l 
(3.6 fim) 


IRAC-2 
(4.5 fxm) 


IRAC-3 
(5.8 /im) 


IRAC-4 
(8.0 fjaa) 


V834 Cen 
EF Eri 
GG Leo 

V347 Pav 
VV Pup 
MR Ser 


6 S4+ - 16 

°- d ^-0.16 

21 + - 05 

u - z± -0.04 

1 cr 7 +0.07 
l -° ' -0.06 

u - d °-0.05 

qfi+ 06 

u - yu -0.06 

o 7-1 +0.09 
°- ' ± -0.09 


5.601°- 

0.5518:8? 

1 97+O.O9 
L - Zl -0.08 

fi^+ 09 

u -°°-0.08 

1 11+O.IO 
± - ±i -0.09 

°- 1 ^-0.11 


r 99+0.13 

n 47+o-os 

u '^' -0.07 

o.93i8;8^ 

o.68l8: 8 

i-oil8:gl 
3.04i8;l? 


2 S9+ - 05 

z - oa -0.05 

71 +0 ' 02 

u - ' i -0.02 
u - zd -0.07 

n 07+0.02 

R7+0- 02 
u -° ' -0.02 

1 so +0 - 03 

±.OU_g 03 


9 oq+0.06 
z -°°-0.06 

n 79+0.02 

u - ' z -0.02 

1 7 +0 - 06 

V - L 1 -0.06 
u -°°-0.02 
u - oz -0.02 

1 44+0-04 

i - 44 -0.04 


2 04+ - 16 

z - u ^-0.16 

66+ 07 

u - uu -0.07 

0.1918:8^ 
o.6ol8:8^ 

n Q4+0.05 

u -°^-0.05 

1 4^+0-12 


2-4818:11 

0.7018:8^ 

1 7 +0 - 03 

v - 1 1 -0.03 
u -°' -0.05 

26 +0 - 04 

u - zu -0.04 

1 17+O.O7 
i - i '-0.07 



Table 3. Polar system parameters from the literature 



Polar 




i 


%D 


Secondary 


M W D 


M 2 


Distance 


B (MG) 




(minutes) 


(°) 


(1000 K) 


Spec. Type 


(M ) 


(M@) 


(pc) 


Dipolc Pole 1, Polc2 



V834 Cen 


101.517 [20] 


50(10) [16] 


14.3(0.9) [1] 


M7.5-8.5 [10] 


0.66(18) [16] 


0.06 [21] 


144j3§ [1] 


23 [4] 




EF Eri 


81.022932(8) [11] 


55(5) [11] 


9.5(0.5) [2] 


M6+/L4-5 [2,7] 


0.6 [2] 


0.055 [9] 


42-120 [2,7] 


13.8 [9] 


16.5, 21 [5] 


GG Leo 


79.87946(7) [3] 


50-60 [3] 


8-12 [3] 


M5+ [3] 


1.13(3) [14] 


0.09 [3] 


100-300 [3] 


20-30 [3] 




V347 Pav 


90.082219(52) [14] 


64-72 [12] 


11.8+°;° [1] 




1.00(6) [14] 


0.16 [13] 


17111 [1] 




15, 20 [12] 


VV Pup 


100.436 [20] 


78(5) [20] 


11.9±°; 6 5 M 


M7 [8] 


0.73(5) [8] 


0.10(2) [8] 


12911? [8] 




32, ■ ■ ■ [18] 


MR Ser 


113.46890(14) [15] 


50(5) [19] 


14.2(0.9) [1] 


M5-8 [17,6] 


0.7 [20] 


0.1 [20] 


irotjg [i] 


28 [20] 


24.5, ■■■ [20] 



= iFerrario et al.l 


1993), [6] 


iMennickent et al.l I 


2004), [111 



lArauic ^Betancore^al ] ll2005l) . [2] 



lHarrison et aU 



lll99ll) . [161 



Beuermann et all ll2000h . [31 = iBurwitz et alj dl998l ). [41 = IFerrario et al.1 lll992h . [5] 
Howell fc Ciardil ll200ll)7[8l = IHowell et alj ll2006rj), [91 = IHowell et alj ll2006ch, [101 = 



. 120051) , [ 71 

=|Piirola et al.1 (Il987l), [121 =|Pptter et al.l ||2000 | ), [131 =|R a msav et a.1.1 Jl996f). [14 
1993al) . [171 = ISchwope et al.l lll993bl ). [181 = IWickramasinghe fc Megeittl Jl98~ 
http : 77cvcat . net | and references therein, [21] = calculated from other parameters (see text) 



iRamsav et al.1 1|2004|)" [15] = | Schwope et al.1 
[191 = IWickramasinghe et al.1 lll99ll) . [201 



Note. — With the exception of distance and secondary star spectral type, the values of these parameters are held constant during our modeling. When 
there is a formal uncertainty or range of values listed here, the fixed value used in the models is the nominal or midpoint parameter value, respectively. 



Table 4. Model Parameters for EF Eri 



Component Parameter a Model 1 Model 2 Model 3 Model 4 Model 5 b Model 6 b 



System: 


d foe) 


105 


105 


105 


105 


105 


105 


SS: 


T2 (K) 


1700 


1700 


1700 


1700 


1790 


1790 




Spec. Type 


L5.0 


L5.0 


L5.0 


L5.0 


L3.5 


L3.5 


CYC: 


B (MG) 




13.8 


13.8 


13.8 




13.8 




??; 




2 


3 


5 




2 




Sum over fields? 




Y 


Y 


N 




Y 




fu,CYc/U,ss (3.6^m) 




1.71 


3.12 


0.34 




1.55 




fu,CYc/U,ss (8.0/im) 




11.01 


7.04 


0.18 




9.34 


CBD: 


Hn (-Rwd) 


73 








73 






T in (K) 


665 








665 






^out (-Rwd) 


7800 








7800 






T out (K) 


20 








20 





a Fixed model parameter values for this target include: -Rwd = 0.0125-R© and R2 = O.IRq. See Table [3] 
for the values of other parameters not listed here. 

b These models refer to the scaled 2MASS data - see text. 



Table 5. Model Parameters for V347 Pav 



Component Parameter 51 Model 1 Model 2 Model 3 Model 4 Model 5 b Model 6 b 



System: 


d (pc) 


171 


171 


171 


171 


133 


133 


SS: 


Spec. Type 


M6.0 


M6.0 


M6.0 


M6.0 


M6.0 


M6.0 


CYC: 


B (MG) 




15,20 


15,20 


15,20 




15,20 




m 




2 


3 


2 




5 




Sum over fields? 




Y 


Y 


N 




Y 




/i/,CYc//i/,SS (3.6/iHl) 




1.06 


0.80 


1.68 




0.62 




fu,CYc/U,ss (8.0/im) 




2.97 


1.48 


0.00 




0.93 


CBD: 


Hn (RWT>) 


148 








148 






T in (K) 


770 








655 






^out (-Rwd) 


19241 








15509 






T out (K) 


20 








20 





a Fixed model parameter values for this target include: -Rwd = O.OO8-R0, R2 = 0.21i? Q , and T2 = 2850 K. 
See Table [3] for the values of other parameters not listed here. 

b These models refer to the scaled 2MASS data - see text. 
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Table 6. Model Parameters for VV Pup 



Component Parameter 3, Model 1 Model 2 



System: d (pc) 100 129 

SS: Spec. Type M7.0 M7.0 

CYC: B (MG) • • • 32 

m ■ ■ ■ 3 

Sum over fields? • • • Y 

fu,CYc/U,SS (3.6/im) ••• 1.13 

fu,CYc/U,SS (8.0/im) ••• 1.35 

CBD: r in (R WD ) 108 

Tm (K) 510 

r ou t (Rwd) 8140 



T out (K) 20 



a Fixed model parameter values for this target include: 
R WD = 0.0105i? o , R 2 = O.15-R , and T 2 = 2670 K. See Table 
[3] for the values of other parameters not listed here. 

Table 7. Model Parameters for V834 Cen 



Component 


Parameter 3, 


Model 1 


Model 2 


Model 3 


System: 


d (pc) 


121 


121 


121 


SS: 


Spec. Type 


M8.0 


M8.0 


M8.0 


CYC: 


B (MG) 




23 


23 




rn 




2 


6 




Sum over fields? 




Y 


Y 




fu,CYc/U,SS (3.6/xm) 




6.57 


6.66 




/^CYC/Z^SS (8.0/xm) 




11.39 


8.37 


CBD: 


n n (-Rwd) 


95 








Tin (K) 


875 








r Q ut (-Rwd) 


14663 








T out (K) 


20 







a Fixed model parameter values for this target include: Rwd = 
0.0115^0, R 2 = O.WRq, and T 2 = 2490 K. See Table [3] for the values of 
other parameters not listed here. 



Table 8. Model Parameters for GG Leo 



Component 


Parameter 3 


Model 1 


Model 2 


kJj O ICJLU. 


d (no) 


900 


900 


SS: 


Spec. Type 


M7.5 


M7.5 


CYC: 


B (MG) 




25 




m 




10 




Sum over fields? 




N 




fv,CYc/fu,SS (3.6//m) 




1.20 




fv,CYc/fu,SS (8.0//m) 




0.00 


CBD: 


r-m (%d) 


185 






T in (K) 


580 






r D ut (%d) 


16504 






T out (K) 


20 





a Fixed model parameter values for this target include: 
% D = O.OO6i? , R 2 = 0.13i?o, and T 2 = 2580 K. See Ta- 
ble [3] for the values of other parameters not listed here. 

Table 9. Model Parameters for MR Ser 



Component 


Parameter 3 


Model 1 


Model 2 


System: 


d (pc) 


134 


134 


SS: 


Spec. Type 


M7.0 


M7.0 


CYC: 


B (MG) 




28 




m 




2 




Sum over fields? 




Y 




fu,CYc/U,SS (3.6/im) 




3.50 




fu,CYc/U,SS (8.0/rni) 




4.79 


CBD: 


r-m (%d) 


122 






Tin (K) 


715 






r D ut (%d) 


14370 






T out (K) 


20 





a Fixed model parameter values for this target include: 
R WD = 0.0LR©, R 2 = O.ISi?©, and T 2 = 2670 K. See Ta- 
ble [3] for the values of other parameters not listed here. 



Table 10. Model Parameters for "Best" Models 



Component Parameter a EF Eri V347 Pav VV Pup V834 Cen GG Leo MR Ser 



System: 


d (dc) 


105 


171 


129 


121 


200 


134 


SS: 


Spec. Type 


L5.0 


M6.0 


M7.0 


M8.0 


M7.5 


M7.0 


CYC: 


B (MG) 


13.8 


15,20 


32 


23 


25 


28 




??; 


3 


2 


3 


6 


5 


2 




Sum over fields? 


Y 


N 


Y 


Y 


Y 


Y 




U,CYc/fu,SS (3.6/zm) 


3.07 


1.20 


1.12 


5.23 


0.94 


3.11 




fu,CYc/U,SS (8.0/um) 


7.06 


0.00 


1.33 


6.56 


1.19 


4.25 


CBD: 


Tin {Rwd) 


73 


148 


108 


95 


185 


122 




T in (K) 


555 


755 


480 


805 


550 


670 




r Q ut (Rwd) 


6129 


18743 


7508 


13119 


15376 


13177 




T out (K) 


20 


20 


20 


20 


20 


20 



Parameters not listed here have the values in the corresponding table from Tables 0HS] and Table [3j 



-45 - 




2468 2468 




2468 2468 




2468 2468 
Wavelength (microns) 



Fig. 1. — Observed spectral energy distributions for our six polars, from 2MASS to IRAC. 












Flux Density (mJy) 



Fig. 2. — Standard deviation of the flux densities in the five BCD frames vs. flux density 
measured in the mosaic for a number of isolated stars in our IRAC fields. Solid line is a best- 
fit polynomial to the data. Dashed line is best fit + rms of residuals (i.e., the la uncertainty 
on the fit). 
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Fig. 3. — Sample model system geometry of EF Eri (face-on, to scale) for the model shown 
in Figure and described in the text. The center of mass (marked with an X) is at the 
origin and the outer Lagrangian points are shown as plus signs. The WD is the small filled 
circle and the secondary star is the large hashed area. The secondary star fills its Roche 
lobe; the Roche lobe of the WD is shown as a dashed line. The solid circle surrounding the 
CV is the inner edge of the circumbinary disk. 
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Fig. 4. — Representative model SEDs for EF Eri showing the observed data (filled circles), 
total model SED (solid line) and model components: the WD (dotted line), the secondary 
star (short dashed line), a circumbinary dust disk (long dashed line), and cyclotron emission 
(dot-dashed line). The plot panels show (a) circumbinary disk model, (b) cyclotron model 
with m = 2 and (c) m = 3, and (d) single- field cyclotron model with m = 5. See text and 
Table H] for details. 
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Fig. 5. — As in Figure HJ but with the 2MASS data for EF Eri scaled by a factor of 1.67 to 
account for potential orbital brightness variations. The plot panels show (a) circumbinary 
disk model and (b) cyclotron model with m = 2. See text and Table |4] for details. 
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Fig. 6. — Representative model SEDs for V347 Pav showing the observed data (filled circles), 
total model SED (solid line) and model components: the WD (dotted line), the secondary 
star (short dashed line), a circumbinary dust disk (long dashed line), and cyclotron emission 
(dot-dashed line). The plot panels show (a) circumbinary disk model, (b) cyclotron model 
with m = 2 and (c) m = 3, and (d) single- field cyclotron model with m = 2. See text and 
Table [5] for details. 
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Fig. 7. — As in Figure EJ but with the 2MASS data for V347 Pav scaled by a factor of 4 to 
account for potential orbital brightness variations. The plot panels show (a) circumbinary 
disk model and (b) cyclotron model with m = 5. See text and Table [5] for details. 
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Fig. 8. — Representative model SEDs for VV Pup showing the observed data (filled circles), 
total model SED (solid line) and model components: WD (dotted line), secondary star (short 
dashed line), circumbinary dust disk (long dashed line), and cyclotron emission (dot-dashed 
line). The plot panels show (a) circumbinary disk model and (b) cyclotron model with 
m — 3. See text and Table [6] for details. 
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Fig. 9. — Representative model SEDs for V834 Cen showing the observed data (filled circles), 
total model SED (solid line) and model components: WD (dotted line), secondary star (short 
dashed line), circumbinary dust disk (long dashed line), and cyclotron emission (dot-dashed 
line). The plot panels show (a) circumbinary disk model and (b) cyclotron model with m = 2 
and (c) m = 6. See text and Table [7] for details. 
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Fig. 10. — Representative model SEDs for GG Leo showing the observed data (filled circles), 
total model SED (solid line) and model components: secondary star (short dashed line), 
circumbinary dust disk (long dashed line), and cyclotron emission (dot-dashed line). At 
maximum, the WD contributes less than 1% of the total observed flux density and is not 
plotted. The plot panels show (a) circumbinary disk model and (b) single-field cyclotron 
model with m = 10. See text and Table M for details. 
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Fig. 11. — Representative model SEDs for MR Ser showing the observed data (filled circles), 
total model SED (solid line) and model components: WD (dotted line), secondary star (short 
dashed line), circumbinary dust disk (long dashed line), and cyclotron emission (dot-dashed 
line). The plot panels show (a) circumbinary disk model and (b) cyclotron model with 
m — 2. See text and Table [9] for details. 
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Fig. 12. — Sample model SEDs that best reproduce the observed data for (a) EF Eri (un- 
sealed), (b) V347 Pav, and (c) VV Pup, showing the observed data (filled circles), total 
model SED (solid line) and model components: WD (dotted line), secondary star (short 
dashed line), circumbinary dust disk (long dashed line), and cyclotron emission (dot-dashed 
line). See text and Table fTUl for details. 



-57- 




Fig. 13. — As in Figure [121 but showing sample "best" models for (a) V834 Cen, (b) GG 
Leo, and (c) MR Ser. See text and Table [10] for details. 



-58- 




Fig. 14. — Speculative model SED for MR Ser showing the observed data (filled circles), 
total model SED (solid line) and model components: WD (dotted line), secondary star (short 
dashed line), and WD circumstellar disk (dot-dot-dashed line - see Appendix |A|) . 



